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1 Introduction 


Since its discovery in 1974, the J meson has provided a useful laboratory for quantita¬ 
tive tests of quantum chromodynamics (QCD) and, in particular, of the interplay of per¬ 
turbative and nonperturbative phenomena. The factorization formalism of nonrelativistic 
QCD (NRQCD) jl] provides a rigorous theoretical framework for the description of heavy- 
quarkonium production and decay. This formalism implies a separation of short-distance 
coefficients, which can be calculated perturbatively as expansions in the strong-coupling 
constant a s , from long-distance matrix elements (MEs), which must be extracted from 
experiment. The relative importance of the latter can be estimated by means of velocity 
scaling rules, i.e. the MEs are predicted to scale with a definite power of the heavy-quark 
(Q) velocity v in the limit v « 1. In this way, the theoretical predictions are organized 
as double expansions in a s and v. A crucial feature of this formalism is that it takes into 
account the complete structure of the QQ Fock space, which is spanned by the states 
n = 25+1 L/y with definite spin S, orbital angular momentum L, total angular momentum 
J, and colour multiplicity c = 1,8. In particular, this formalism predicts the existence of 
colour-octet (CO) processes in nature. This means that QQ pairs are produced at short 
distances in CO states and subsequently evolve into physical, colour-singlet (CS) quarko- 
nia by the nonperturbative emission of soft gluons. In the limit v —» 0, the traditional CS 
model (CSM) |2j.[3j is recovered. The greatest triumph of this formalism was that it was 
able to correctly describe [f||| the cross section of inclusive charmonium hadroproduction 
measured in pp collisions at the Fermilab Tevatron [E|, which had turned out to be more 
than one order of magnitude in excess of the theoretical prediction based on the CSM. 

In order to convincingly establish the phenomenological significance of the CO pro¬ 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and uN deep-inelastic 
scattering (DIS), e + e _ annihilation, 77 collisions, and b -hadron decays may be found in 
the literature; see Ref. |7[] and references cited therein. Furthermore, the polarization 
of charmonium, which also provides a sensitive probe of CO processes, was investigated 
[§|J||TCJ. None of these studies was able to prove or disprove the NRQCD factorization 
hypothesis. 

In this paper, we revisit J/iJ) inclusive production in ep DIS. In order to avoid kinematic 
overlap with diffractive production, which ca nn ot yet be reliably described within purely 
perturbative QCD, we require that the J/ij) meson be produced in association with a 
hadron jet j, i.e. we consider the process e + p —> e + J/ijj + j + A", where A" denotes the 
proton remnant. In this way, the inelasticity variable z, which measures the fraction of 
the virtual-photon ( 7 *) energy transferred to the J/rp meson in the proton rest frame, can 
take values below unity, away from the endpoint z = 1 , where diffractive production takes 
place. At the same time, the J meson can acquire finite transverse momentum p ^ in 
the 7 *p centre-of-mass (CM) frame, and the hadronic system X' consisting of the jet j and 
the proton remnant X can acquire finite mass M x >- By the same token, diffractive events 
can be eliminated from the experimental data sample by applying appropriate acceptance 
cuts on z, Pt, or M x >. Another possibility to suppress the diffractive background at z 1 
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would be to require that the photon virtuality Q 2 be sufficiently large [|1T|] . However, then 
also the bulk of the nondiffractive signal would be sacrificed. 

The leading CS ME of the J/-0 meson is fS'jits leading CO ones are 

< 'O* [tfj 0 ]), ( O* [ 3 5f } ]), and (O^ [ 3 Pj 8) ]), with J = 0,1,2. At LO, we are thus led 
to consider the partonic subprocesses e + a —> e + cc[n] + a, where a = g,q,q and 
n = 3 S? , 3 S[ 8 \ 3 Pf. Here, q runs over the light-quark flavours u, d, and s. Notice 


that e + q(q ) —> e + cc 


^ +q(q) is forbidden because the charm-quark line is connected 
to one gluon, which transfers colour to the cc pair. Representative Feynman diagrams are 
depicted in Fig. [Tj. The corresponding cross sections are of 0(a 2 a 2 ). 

There are several motivations for our study oie + p—te+J/ip + j + X in DIS. For 
one thing, the HI Collaboration measured various distributions of this process at DESY 
, P~4| , which wait to be confronted with appropriate theoretical predictions. 


HERA 

This allows for a particularly clean test of the NRQCD factorization hypothesis, since 
the large photon virtuality Q 2 ensures that perturbative QCD is applicable and that the 
resolved-photon contribution, which suffers from our imperfect knowledge of the parton 
density functions (PDFs) of the photon, is greatly suppressed. Furthermore, at least 
within the CSM, this process provides a good handle on the gluon PDF of the proton, 
which is less precisely known than the quark ones. In NRQCD, however, the extraction 
of the gluon PDF is somewhat aggravated by the presence of the CO channels with a 
quark or antiquark in the initial state and by the uncertainties associated with the CO 
MEs. Fortunately, detailed inspection reveals that the relative importance of the quark- 
and antiquark-induced channels is greatly damped in the HERA regime. We return to 
this point at the end of Section |3|. 

On the other hand, in the case of J/-0 inelastic photoproduction, with Q 2 « 0, NRQCD 
with CO MEs tuned [|] to fit the Tevatron data |J predicts |I5| , [T6[ at leading order (LO) 
a distinct rise in cross section as z —> 1, which is not observed by the HI Jl7j and ZEUS 
m Collaborations at HERA. This CO charmonium anomaly has cast doubts on the 
validity of the NRQCD factorization hypothesis, which seems so indispensible to inter¬ 
pret the Tevatron data in a meaningful way. Although there are several interesting and 
promising ideas how to reconcile this data with the NRQCD prediction, e.g. by including 
dominant higher-order effects [IP], by introducing nonperturbative shape functions that 
resum higher-order corrections related to the kinematics of soft-gluon radiation and to 
the difference between the partonic and hadronic phase spaces [^0[, or by endowing the 
partons inside the proton with intrinsic transverse momentum (hr) a it is of great in¬ 
terest to fold out if this anomaly persists, at LO and without resorting to shape functions 
or k T effects, if Q 2 is increased to large values. 

On the theoretical side, the cross sections of the partonic subprocesses under considera¬ 
tion here constitute an essential ingredient for the calculation of the next-to-leading-order 
(NLO) corrections to the inclusive cross section of the reaction e + p—►e + J/^ + A’ in 
DIS, which, at LO, proceeds through the partonic subprocesses e + g —> e + cc[n\ with 
n = 1 4 8) , 3 ^j 8) - The cross sections of the latter are of 0(a 2 a s ) and may be found in 
Eq. (3) of Ref. |22|. In fact, integrating over the phase space of the massless final-state 
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parton a, one obtains the real radiative corrections. If a = g and n = ■'S'®, 3 Pj 8 \ then 
they suffer from both infrared (IR) singularities and collinear ones associated with the in¬ 
coming gluon. The latter are factorized, at some factorization scale M, and absorbed into 
the bare gluon PDF of the proton, so as to render it renormalized. The IR singularities 
cancel when the real radiative corrections are combined with the virtual ones. Finally, 
the ultraviolet (UV) radiative corrections contained in the latter are removed by renor¬ 
malizing the couplings, masses, wave-functions, and non-perturbative MEs appearing in 
the LO cross section of e + p —> e + J/if} + X. 

Finally, the literature contains mutually inconsistent formulas for the cross section of 


+ gm 


27]. On the other hand, 


the CS partonic subprocess e + g —> e + cc 
there is only one paper specifying analytic results for the cross sections of the CO partonic 
subprocesses enumerated above [|27]], so that an independent check seems to be in order. 
We anticipate that we disagree with all published CS and CO formulas, except with the 
one referring to e + g — > e + cc + g |27[ . 

This paper is organized as follows. In Section |2], we present, in analytic form, the 
cross sections of the partonic subprocesses e + a —»■ e + cc[n} + a enumerated above 
to LO in NRQCD and explain how to calculate from them the total cross section of 
e + p —> e + J/'if) + j + X in DIS and several distributions of phenomenological interest. 
Lengthy expressions are relegated to the Appendix. In Section [3], we present our numerical 
results and compare them with recent HI data [|l^]. Our conclusions are summarized in 
Section f|. 


2 Analytic results 


In this section, we present our analytic results for the cross section of e+p —■> e+J/ip+j+X 
in DIS. We work at LO in the parton model of QCD with rif = 3 active quark flavours 
and employ the NRQCD factorization formalism [I| to describe the formation of the J/rj} 
meson. We start by defining the kinematics. As indicated in Fig. |2j. we denote the four- 
momenta of the incoming lepton and proton and the outgoing lepton, J/ip meson, and 
jet by k, P, k', p^,, and p', respectively. The parton struck by the virtual photon carries 
four-momentum p = xP. We neglect the masses of the proton, lepton, and light quarks, 
call the one of the J/ij) meson M,/,, and take the charm-quark mass to be m c = M^/2. In 
our approximation, the proton remnant X has zero invariant mass, Mj- = (P — p) 2 = 0. 
The CM energy square of the ep collision is S = (k + P) 2 . The virtual photon has four- 
momentum q = k — k ', and it is customary [2S| to define Q 2 = —q 2 > 0 and y = q-P/k-P, 
which measures the relative lepton energy loss in the proton rest frame. The inelasticity 
variable, which was already mentioned in the Introduction, is defined as z = p^ ■ P/q ■ P. 
The system X' consisting of the jet j and the proton remnant X has invariant mass 
square M\, = (q + P — p^) 2 = (1 — x)y( 1 — z)S. Other frequently employed variables 
p8| are Bjorken’s variable xb = Q 2 /{ 2 g ■ P) = Q 2 /(yS ) and the CM energy square 
W 2 — (q + P ) 2 = yS — Q 2 . As usual, we define the partonic Mandelstam variables as 
s = (q +p) 2 = xyS - Q 2 , t = (q -p </,) 2 = —xy{ 1 - z)S , and u = (p - p^) 2 = M 2 - xyzS. 
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By four-momentum conservation, we have s + t + u = M 2 — Q 2 . In the CM frame, 
the J/'ip meson has transverse momentum and rapidity 


\jt (su + Q 2 

Pt = 


s + Q 2 


o ln “ 

^ s 


(Ml - u) 


(Ml - i) + QHll 


1 , w 2 

H— In-. 

2 2 s' 


( 1 ) 

( 2 ) 


respectively. Here and in the following, we denote the quantities referring to the j*p CM 
frame by an asterisk. The second term on the right-hand side of Eq. (|2|) originates from 
the Lorentz boost from the y*a CM frame to the one. Here, y ^ is taken to be positive 
in the direction of the three-momentum of the virtual photon, in accordance with the 
convention of Refs. 


12,13,14 


In the parton model, the proton is characterized by its PDFs f a / P (x, M), and, at LO, 
an outgoing parton may be identified with a jet. Thus, we have 


da(e + p^>e + J/ip + j + X ) — J dx'y fa/ P (x, M)da(e + a —> e + J /-0 + a), 


( 3 ) 


where a = g,u, w, d , d, s, s. Furthermore, according to the NRQCD factorization formal¬ 
ism, we have 

dcr(e + a —> e + J /-0 + a) = y dcr(e + a —> e + cc[n] + a), (4) 


where, to LO in v, n = 3 S[ 8 \ 3 Pj 8 \ 

Decomposing the transition-matrix element of the partonic subprocess e + a — > e + 
cc[n\ + a into a leptonic part, T M (e —> e + 7 *) = — (e / q 2 )u(k , )'y fl u(k), and a hadronic one, 
T M (y* + a —> cc[n] + a), from which the virtual-photon leg is amputated, we can write its 
cross section as 

lie 2 

da(e + a -*• e + cc[n] + a) = AV)#^PS 3 (fc + p; k',p^,p'), (5) 

where 1V 9 = (IV 2 — 1) and iV g = Ng = N c = 3 are the colour multiplicities of the partons a, 
e is the electron charge magnitude, and the hadronic tensor H^ is obtained by summing 
the absolute square of T M (y* + a —> cc[n] + a) over the spin and colour states of the 
incoming and outgoing partons a. Here and in the following, we employ the Lorentz- 
invariant phase-space measure 


dPS n (p;pi, ...,p n ) 


(27t) 4 £ (4) [p-ypi ) n 


2=1 


2=1 


d 3 pi 

(2vr) 3 2pC 


( 6 ) 


The first factor in Eq. (|5|) stems from the flux and the second one from the average over 
the spin and colour states of the incoming particles. Integrating over the azimuthal angle 
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of the outgoing lepton, we may simplify Eq. (|5|) to become 

11 a 


dc r(e + a —> e + cc[n] + a) = 


2xS' 41V„ 27r 


dy dQ 2 

V~Q 2 


—L^H^ ^dPS 2 (q + p-p^p'), (7) 


where a = e 2 /( 47 r) is Sommerfeld’s fine-structure constant and |29 


jjiv = i + (i -y? c w? _ 4 ( x -y) 


y 


y 


A1V 

e L i 


( 8 ) 


with 


+ 


1 («y+pV)- 7 ^3jp , y', 


£f = 7 9 


9 -V 

q-p 

_fiiy _ 


(q-pf 


-p 


Q 


Q 


-p 


q-p 
cV-v _ n ^ — 


( 9 ) 


is the leptonic tensor. We have = q^e^ = 0, ^ = —2, and e'/ = —1. Furthermore, 


= eg + ef = + 


g /J g' y 


( 10 ) 


is the polarization tensor of an unpolarized spin-one boson with mass q 2 . In the y*a 
CM frame, where q M = (g°, 0, 0, q 3 ) and p M = (g 3 , 0, 0, — g 3 ), we have = diag(0,1,1, 0) 
and dj" = (1 /g 2 )(g 3 , 0 , 0 , g°) / 9 (g 3 , 0 , 0 , q°) u , so that e? and refer to transverse and 
longitudinal polarization, as indicated by their subscripts. Since the hadronic current is 
conserved in QED, we have q^H^ = 0, which leads to a further simplification as Eq. (|S]) 
is inserted in Eq. (^). 

It is interesting to study the photoproduction limit, by taking Q 2 —> 0 in Eq. (0). 
This provides us with a powerful check for our results by relating them to well-known 
results in the literature The differential cross section of the partonic process 

7 + a —> cc[n] + a reads 


dcr (7 + a —> cc[n] + a) = — 


-g^)H^\ Q 2 =0 dPS 2 (q + p-,PrpiP')- 


2 s 4 N a 

Comparing Eqs. (0) and (0), we thus obtain the master formula 

a 1 + (1 -y) 2 


( 11 ) 


Q 2 d 2 a 


e + cc[uj + a) = — 

Z7T 


cr (7 + a —> cc[n\ + a). ( 12 ) 


A similar relationship between the DIS process e + p —> e + J/i/j + X and the photopro¬ 
duction one 7 + p —> J/ip + X may be found in Eq. (4) of Ref. |22| . 

We evaluate the cross sections of the relevant partonic subprocesses e+a —*• e+cc[n]+a 
from Eq. ([?]) applying the covariant-projector method of Ref. |3(J. Our results can be 
written in the form 


d 3 a / _ r n \ a r , 

- -(e + a —> e + cc\n\ + a) = — r a \n\ 

dydQ 2 dV 1 J ; 27t 


1 + 


yQ 2 


y 


, ( 13 ) 
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where F a [n ], T a [n], and L a [n] are functions of s, t , u, and Q 2 , which are listed in the 
Appendix. They are finite for Q 2 = 0. We combined the results for n = 3 Pj 8 \ with 
J = 0,1, 2, exploiting the multiplicity relation 



(2 J + 1) (O* [ 3 P 0 (8) ] ) 


(14) 


which follows to LO in v from heavy-quark spin symmetry. We recover the well-known 
cross sections of the corresponding CS [|J and CO |||I7j processes of photoproduction by 
inserting Eq. fll3D in Eq. (|l^) , as 


^(7 + a cc[n] + a) = F a [n]T a [n] | Q 2 =0 . 

At this point, we should compare our analytic results with t 
for the cross section of the CS partonic subprocess e + g —> e + cc 
in Refs. 


24 


(15) 

re literature. Formulas 
w + g may be found 
In Ref. |27[], also cross section formulas for the CO partonic 


e + cc[n] + a with a = g,q,q and n = iS®, 3 S[ 8 \ 3 Pj 8 ' > are listed. We 

15 .( 8 ) 


subprocesses e + a 

agree with the result for e + g —> e + cc 
results. I 11 particular, the results of Refs, 
reproduce the well-established formulas of Refs. p]M[l 6 | in the photoproduction limit. We 


+ g |27|], but we disagree with all the other 
and the residual results of Ref. [^7| fail to 


also remark that Eqs. (A 21 ) and (A35)-(A37) of Ref. p7[ suffer from mass-dimensional 


inconsistencies. Furthermore, we only find agreement with the result in Eq. (4) of Ref. [^4 
if we flip the overall sign of L a [n\ in Eq. (jh|). However, this causes the cross section in 
Eq. (jT3|) to turn negative in certain regions of phase space, e.g. at y = 0.5, Q 2 = 25 GeV 2 , 
s = 100 GeV 2 , and t — —10 GeV 2 . Similarly, we only agree with the result in Eq. (2.49) 
of Ref. [Q if we include an overall factor of Ana on the right-hand side of Eq. ([T5|) and 
halve L a [n\. As for the cross sections of the CO partonic subprocesses e + g —► e + cc[n] 
with n = ^Sq 8 ) , 3 Pj 8 \ which do not enter our analysis, we agree with Eq. (3) of Ref. [22 . 

The cross sections of e + g —> e + cc[n] + g with n = 3 Pj 8 ^ exhibit collinear 

singularities in the limit t —*■ 0. According to the factorization theorem of the QCD- 
improved parton model, the limiting expressions must coincide with the e + g —» e + cc[n] 
cross sections multiplied by the spacclike g —> g splitting functions. This provides 


another nontrivial check for our results, and, among other things, this fixes the overall 
factor of L a [n\. 

Inserting Eq. (|4j) in Eq. (|3j) and including the maximum boundaries of the integrations 
over x and t, we obtain 


d 2 


a 


dy dQ 2 


x 


(e + p—^e + J/'i/j + j + X) = 

d 3 a 


r 1 r 0 

dx I 

12 1 s,r2 -i 


dt 




dy dQ' 2 dt 


(Q 2 +M|)/(yS) J-(s+q2)(s-MQ/s 

(e + a —> e + cc[n] + a), (16) 


where (cl 3 a/dy dQ 2 dt) (e + a —> e + cc[n] + a) is given by Eq. (|T3|). The kinematically 
allowed ranges of y and Q 2 are M 2 /S < y < 1 and 0 < Q 2 < yS — M 2 , respectively. In 
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order to avoid the collinear singularities mentioned above, we need to reduce the upper 
boundary of the t integration in Eq. ([16]) . This cut-off should, of course, depend on 
variables that can be controlled experimentally. It is convenient to introduce an upper 
cut-off, below unity, on the inelasticity variable z. As explained in the Introduction, 
such a cut-off also suppresses the diffractive background. The distributions in y and 
Q 2 can be evaluated from Eq. (|16|) as it stands. It is also straightforward to obtain 
the distributions in z , Xb, W, p?, and y^, by accordingly redefining and reordering the 
integration variables in Eq. (|l^). The distribution in the J/p azimuthal angle </>* in the 
7 *p CM frame is constant. 

The evaluation of the distributions in the J/p transverse momentum p?, rapidity 
y^, and azimuthal angle (j) in the HERA laboratory frame is somewhat more involved. 
Choosing a suitable coordinate system in the CM frame, we have 


= 


S-Q 2 

2W 


( 1 ) 

- 

( 1 \ 

sin p* 

0 

(k »Y - S ~ W 
’ ^ ’ 2 IT 

sin 9* 

0 

\ COS if)* / 


\ cos 6* J 


(gY 


(W 2 — Q 2 \ 
1 0 

2 W 0 

\W 2 + Q 2 ) 


(p*y = 


W 2 + Q 2 
2 W 


( 1 \ 


0 

0 ’ 
V-i/ 



/ mtp cosh y ^ \ 
Pt cos (j )* 

Pt sin </>* 

\ nip sinh y£ / 


( 17 ) 


where cos^» = 2SM-'7[( S - Q 2 )(W 2 + Q 2 )] - 1 , cos9* = 1 - 2SQ 2 /[(S- W 2 )(W 2 + Q 2 )], 
and M 2 + (pp) 2 is the J/p transverse mass. On the other hand, in the laboratory 

frame, we have 



f 1 ) 


( 1 ^ 

k» = E e 

0 

0 

, (k'y = K 

sin 6 

0 


\l) 


V cos 9 J 


( q ° ) 


( 1 \ 

—q sin p 

, = E v 

0 

0 

1 P 

0 

\ geos p / 


\-l) 


/ m T cosh y,e \ 
Pt cos (j) 

Pt sin (j) 

V m T sinh y^ ) 


( 18 ) 



























where E e and E p are the lepton and proton energies, respectively, E' e = [(S' — W 2 — 
Q 2 )/E p + Q 2 /E e \/ 4, cos 6 = 1 - Q 2 /(2E e E ' e ), g° = [(W 2 + Q 2 )/ E p - Q 2 / E e ]/4, q = 
\JQ 2 + (g 0 ) 2 , cos ip = [{W 2 + Q 2 )/E p + Q 2 /£ e ]/(4g), and m T = is the J/-0 

transverse mass. Notice that y^ is taken to be positive in the direction of the three- 
momentum of the incoming lepton. Without loss of generality, we may require that 
0 < 0*, 0 < 7T, for, otherwise, we can achieve this by rotating the respective coordinate 
systems by 180° around the z axis. We can then evaluate px, y^, and 0 from ptp, y^, and 
0 * as 


Pt 

y^p 

COS 0 


\J (; Pt) 2 + A (A - 2p* T cos 0*), 


^ VSWrriT 

p* T cos 0* — A 

Pt 


+ bn —, 

2 E p 


(19) 

( 20 ) 

( 21 ) 


where 


A = 


m? exp (yA) sin-0 


lQ 2 (S -W 2 -Q 2 ) 


mjexp (yj). 


( 22 ) 


1 + cos ip * V 5W 2 
The third term on the right-hand side of Eq. (|20|) stems from the Lorentz boost from the 
ep CM frame to the laboratory one. Since px, yp, and 0 depend on 0*, the integration over 
0* in Eq. ([16]) is no longer trivial, and we need to insert the symbolic factor (1/27 t) / 0 27r d(f>* 
on the right-hand side of that equation. 

For future applications, we also present compact formulas that allow us to determine 
Px, y and 0*, once px, y-ip, and 0 are given. In fact, Eqs. (|T9| ) and (^1]) can be straight¬ 
forwardly inverted by observing that the quantity A defined in Eq. (|22|) can be expressed 
in terms of mx and yp by substituting Eq. (|2U|), the result being 


A 


rn T exp(y,/,) 
W 2 + Q 2 



Q 2 (S -W 2 -Q 2 ). 


(23) 


Having obtained p^, we can then evaluate y^ from Eq. (^0|). For the reader’s convenience, 
we collect the relevant formulas here: 


Px 

4 

COS 0* 


x + A{A + 2 p T cos 0 ), 

, VSWm T 1, E v 

\yi - 1 — m — 

{W 2 + Q 2 )m* T 2 E e ' 
cos 0 + A 


(24) 


3 Numerical results 

We are now in a position to present our numerical results. We first describe our theoretical 
input and the kinematic conditions. We use m c = (1.5 ± 0.1) GeV, a = 1/137.036, and 
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the LO formula for a s nf \n ) with rif = 3 active quark flavours [ESI. As for the proton 
PDFs, we employ the LO set by Martin, Roberts, Stirling, and Thorne (MRST98LO) 
0 , with asymptotic scale parameter A^ 4 ) = 174 MeV, as our default and the LO set 
by the CTEQ Collaboration (CTEQ5L) [[Rj, with A^ = 192 MeV, for comparison. The 
corresponding values of A^ 3 ^ are 204 MeV and 224 MeV, respectively. We choose the 
renormalization and factorization scales to be n = M = £ \JQ 2 + M 2 and vary the scale 
parameter £ between 1/2 and 2 abound the default value 1. We adopt the NRQCD 
MEs from Table 1 of Ref. [TOj. Specifically, they read ^ = (1.3 ± 0.1) GeV 3 , 

(p* [ 3 S{ 8) ]) = (4.4 ± 0.7) x 10 “ 3 GeV 3 , and M% A = (8.7 ± 0.9) x 10 " 2 GeV 3 for set 
MRST98LO, where 

" ' + A (O* fef)) . (25) 


M* = (O* (‘S, 


( 8 ) 

0 


mi 


The corresponding values for set CTEQ5L are (1.4 ± 0.1) GeV 3 , (3.9 ± 0.7) x 10 3 GeV 3 , 
and ( 6.6 ± 0.7) x 10 -2 GeV 3 , respectively. Since Eq. (|T^) is sensitive to a different lin¬ 
ear combination of (O^ (^Sq 8 ^ and (O^ ( 3 To 8 ^)) than appears in Eq. (0), we write 

( O* ( 1 S'J 8) )} = nMf and (o* ( 3 P 0 (8) )} = (1 - k) (m 2 /r) Mf and vary k between 0 and 
1 around the default value 1 / 2 . In order to estimate the theoretical uncertainties in our 
predictions, we vary the unphysical parameters £ and n as indicated above, take into 
account the experimental errors on m c and the default MEs, and switch from our default 
PDF set to the CTEQ5L one, properly adjusting A® and the MEs. We then combine 
the individual shifts in quadrature. 

The HI data on J/-0 inclusive and inelastic production in DIS [|1^JT^| were taken in 
collisions of positrons with E e = 27.5 GeV and protons with E p = 820 GeV in the HERA 
laboratory frame, so that a fS = 2 E e E v = 300 GeV, and they refer to the kinematic 
region defined by 2 < Q 2 < 80 GeV 2 , 40 < W < 180 GeV, and z > 0.2. In Ref. fTBfl , the 
acceptance cut M x > > 10 GeV on the invariant mass of the hadronic system X' produced 
in association with the J /-0 meson was imposed in order to exclude the contribution of 
J/-0 elastic production, which, to a large extent, is due to diffractive processes. Notice 
that this cut allows for z to be as large as z max = 1 — M\,( w 2 ,^ — M 2 ^ & 0.997 


and for 


to be as small as 


= (1 


(Qr 


+ M 2 


0.040 GeV 2 . A more 


conservative way to eliminate the domain of J/^ elastic production is to directly impose 


the cut z < 0.9 |12|, which ensures that 
time, a sufficiently large lower bound on 


is in excess of 


1.3 GeV 2 . At the same 


is requisite in order to screen Eq. ( |T7)| ) from the 

collinear singularities in the and 3 Pj 8 ^ channels, mentioned in Section |2j. and, thus, 
to keep our theoretical predictions perturbatively stable. In the following, we, therefore, 
use the HI data with the cut z < 0.9, as presented in Ref. [O, for comparisons. 


In Figs. 


we confront the measured Q 2 


Pti 


z, y *, and W distributions, respec¬ 


tively, with our NRQCD predictions. For comparison, we also show the corresponding 
GSM predictions. In each case, the theoretical errors, evaluated as explained above, are 
indicated by the hatched areas. Instead of presenting our theoretical predictions as con- 
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tinuous curves, we adopt the binning pattern encoded in the experimental data, so as to 
facilitate quantitative comparisons. As for the Q 2 and Px distributions, the experimental 
data agrees rather well with the NRQCD predictions, both in normalization and shape, 
while it significantly overshoots the CSM predictions. In the case of the z distribution, the 
NRQCD prediction, in general, agrees better with the experimental data than the CSM 
one as far as the normalization is concerned. As for the shape, however, the experimental 
measurement favours the CSM prediction, while the NRQCD one exhibits an excess at 
large values of z, which is familiar from J/ij) inclusive photoproduction |T5|JT6fl . As in 


the latter case, this rise in z is chiefly due to the and 3 Pj 8) channels. In the case 
of the y* distribution, the NRQCD prediction nicely agrees with the experimental data 
for y* < 3, while it appreciably overshoots the latter in the very forward direction. The 
CSM prediction roughly agrees with the experimental y * distribution at the endpoints, 
while it significantly falls short of the latter in the central region. In the case of the W 
distribution, the experimental data mostly lie in the middle between the NRQCD and 
CSM predictions, slightly favouring the former. We conclude from Figs. [|-[7]that the HI 
data III tends to support the NRQCD predictions, while, in general, it overshoots the 
CSM predictions. In Figs. |S] and |S| we present our NRQCD and CSM predictions for 
the (Pt) 2 and y distributions, respectively, although there are no experimental data to 
compare them with. 

At this point, we should compare our numerical results for the cross section of e + p —» 


e + J/ip + j + X with the ones presented in Table II of Ref. 
adopt the theoretical input and kinematic conditions from Ref. 


221 . 


22 


To this end, we 
Specifically, the 


p2| evaluated ai nf \p) with rif — 4 and A^ = 130 MeV, employed 

|k| , took the NRQCD MEs to be 
= 1.12 xl0~ 2 GeV 3 , 


authors of Ref. 

the LO proton PDF set by Gliick, Reya, and Vogt 
(O^ [ 3 R( 1} ]) = 1.1 GeV 3 , ( O* [^ 8) ]) = 1 x 10“ 2 GeV 3 , (p^ [ 3 S{ 8) 

and 3 P 0 (8) ^ /m 2 = 5 x 10~ 3 GeV 3 , and required that 30 < W < 150 GeV. All their 
other choices, except for the cuts on Q 2 , px, Px, and z, coincide with our default settings. 
The outcome of this comparison is presented in Table [[]. We are unable to determine the 
source of discrepancy. 


Table 1: Comparison of our results for the CS and CO contributions to cr(e + p 
e + J/ + j + X) in DIS with the ones of Ref. 


22 


Type 

Cuts 

Ref. H 

Our result 

CS 

CS 

CS 

CO 

Q 2 > 4 GeV 2 

Q 2 ,p 2 x> 4 GeV 2 

Q 2 > 4 GeV 2 , G^) 2 > 2 GeV 2 , z < 0.8 
Q 2 > 4 GeV 2 , G4) 2 > 2 GeV 2 , z < 0.8 

89 pb 
40 pb 
13 pb 

8 pb 

107 pb 
62 pb 
24 pb 
16 pb 


Before the advent of the NRQCD factorization formalism |l]], one of the major moti- 
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vations to study J/ip inclusive production in ep DIS was to extract the gluon PDF of the 
proton fg/ p (x, M ) )[24|.[2^|1 . In fact, to LO the CSM, the only contributing partonic subpro¬ 


cess is e + g —> e + cc 


[ S si 


(i) 


is well determined from the partial width 


+ g, and [ 3 S( 1} " 

of the J/ip decay to lepton pairs. Furthermore, x is experimentally accessible through 
the relation x — 1 — M \,/ [y (1 — z)S], and M = ^jQ 2 + is a plausible choice. To LO 

in NRQCD, this task is somewhat impeded by the presence of CO partonic subprocesses 
with quarks or antiquarks in the initial state and by the presently still considerable uncer¬ 
tainties in the CO MEs. We observe from Figs. M that the CO contributions lead to a 
dramatic increase in cross section relative to the CSM predictions. In fact, if we consider 
the total cross section evaluated with the HI acceptance cuts 0- then the CSM contri¬ 
bution only makes up 25% of the full NRQCD result. Thus, the present uncertainties in 
the CO MEs do constitute a serious problem. However, the quark- and antiquark-induced 
CO subprocesses only yield a minor contribution, less than 5% of the total cross section. 
Thus, J/ip inclusive production in ep DIS remains to be a sensitive probe of the gluon 
PDF of the proton if we pass from the CSM to NRQCD. 


4 Conclusions 


We provided, in analytic form, the cross sections of the partonic subprocesses e + a —> 
e + cc[n] + a, where a = g,q,q and n = ^j 1 *, U;®, 3 <Si , 3 Pf \ to LO in the NRQCD 
factorization formalism |l|. Using these results, we then studied the cross section of 
e + p —» e + J/ip + j + X in DIS under HERA kinematic conditions and compared 
various of its distributions with recent HI data 
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from which the contribution of elastic 
J/ip production was separated by a suitable acceptance cut on 0 . It turned out that, 
in general, the experimental data agrees reasonably well with our NRQCD predictions, 
while they tend to disfavour the CSM ones. However, a familiar feature of J/ip inclusive 
photoproduction in ep scattering at HERA [TROR] carries over to DIS, namely that the 
cross section predicted to LO in NRQCD exhibits a distinct rise as z —» 1, which is absent 
in the experimental z distribution. 

At this point, it is still premature to jump to conclusions concerning the experimental 
verification or falsification of the NRQCD factorization hypothesis. On the one hand, 
the theoretical predictions for J/ip inclusive production in pp scattering at the Tevatron 
[40] and in ep DIS at HERA are of LO in a s and v, and they suffer from considerable 
uncertainties, mostly from the scale dependences and from the lack of information on the 
CO MEs. On the other hand, the experimental errors are still rather sizeable P JT4JT3| . 
The latter will be dramatically reduced with the upgrades of HERA and the Tevatron, and 
with the advent of CERN LHC and hopefully a future e + e _ linear collider such as DESY 
TESLA. On the theoretical side, it is necessary to calculate the NLO corrections to the 
hard-scattering cross sections and to include the effective operators which are suppressed 
by higher powers of v. 
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A Partonic cross sections 


In this Appendix, we present analytic expressions for the coefficients F a [n], T a [n ], and 
L a [n] appearing in Eq. (|T3|). In order to compactify the expressions, it is useful to intro¬ 
duce the Lorentz invariants s = 2q ■ p, t = —2 p ■ p ', and u = —2 q ■ p\ which are related to 
the partonic Mandelstam variables by s = s + Q 2 , t = t, and u = u + Q 2 , respectively. In 
the following, e q is the fractional electric charge of quark q. 


e + q(q) 
e + q(q) 


cc 


cc 


+g(g): 


F = T = L = 0. 


(A.l) 


1 ^( 8 ) 


+ q(q ) : 


F = 


— 167 T 2 e 2 aa 2 


3 M^sH(s + u ) 2 ’ 

T = 2 QH 2 + 2 Q 2 st(s + u) + s 2 (s 2 + u 2 ), 
L = —2t(Q 2 t + su ). 


(A.2) 


+ <?(<?) —► cc 3 5{ 8) + q(q): 


F = 


—8n 2 e 2 aa 2 


9 M 3 s 4 (Q 2 -s) 2 (Q 2 -u) 2 ’ 

T = 2Q 6 t 2 (2s + t) + 2 Q 4 s[s 2 u — st(3t — 2 u) — 2t 3 ] + Q 2 s 2 [s 2 (t — 2 u) 

+ 2 s(t 2 — Atu — u 2 ) + t(2t 2 — 2 tu — u 2 )\ + s 3 u(s 2 + 2 st + 21 2 + 2 tu + u 2 ), 

L = — 2{Q 2 — s) 2 t(s + t) 2 . (A.3) 


q{q) 


cc 


3p m 


+ q(q): 


F = 


647r 2 e^o;Q;g 

3 + u) 
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T = 8Q d t[2s 2 + st + t(2t + m)] — 2Q 4 [4s 4 + 12 s 3 t + s 2 (19t 2 + 8 tu + 4m 2 ) + 2st(6t 2 
+ tu — 2 u 2 ) + t 2 (8t 2 + 12 tu + 7m 2 )] + 2Q 2 s[6s 4 + s 3 (7t + 6m) + s 2 (4t 2 + 3 tu 
+ 6m 2 ) — su(8t 2 + 3 tu — 6m 2 ) — tu(8t 2 + 12£m + 7m 2 )] — s 2 (s + m)[7s 3 + s 2 (12t 
+ 7 m) + s(8t 2 + 16£m + 7m 2 ) + u(8t 2 + 12 tu + 7m 2 )], 

L = 8 QH[s 2 — st — t(2t + m)] — 2Q 2 [2s 4 + 4 s 3 t + s 2 (5t 2 + 8£m + 2m 2 ) — 2st(6t 2 
+ tu — 2m 2 ) — t 2 (8t 2 + 12tM + 7 m 2 )] + 2 s(s + m)[2s 3 + 2s 2 t + s(8t 2 + 5£m 
+ 2m 2 ) + t(8t 2 + 12tM + 7 m 2 )] . (A.4) 

e + g —> cc 3 5i^ + g: 

647r 2 e 2 aa 2 

27 M^s 4 (s + t) 2 (s + u) 2 (t + m ) 2 ’ 

T = —4Q 6 t 2 (s 2 + t 2 ) + 2Q 4 t[s 3 (3t — 2m) + 3s 2 t(t + m) + 2st 2 (t — m) + 2t 3 (t + m)] 

— 2Q 2 s[s 3 (t — m) 2 — 2s 2 tu{t + m) — st 2 u{2t — u) — 2 t 3 u(t + m)] + 2s 2 [s 3 (t 2 
+ tu + M 2 ) + S 2 (t + m) 3 + stu{t 2 + 3£m + M 2 ) + t 2 U 2 (t + m)], 

L = —2Q 4 t 2 (s 2 — 2 t 2 ) + 2 Q 2 t{s 2 — 2 t 2 )[s(t — u) + t{t + m)] — s[s 3 (t 2 + m 2 ) 

+ 2 s 2 t 2 {t + m) + st 2 {t 2 + 6tu + m 2 ) + 4t 3 M(t + m)]. (A.5) 

e + 9 -> cc pSj 8) ] + g: 

_ 247r 2 e 2 aa 2 

M^sH(s + t) 2 (s + m) 2 (£ + m ) 2 ’ 

T = 2 Q 6 t 3 u 2 + 2Q 4 st 2 u[s(t + m) + t 2 + tu + 2m 2 ] + Q 2 s 2 t[s 4 + 2s 3 (t + m) 

+ 3s 2 (t + m ) 2 + 2s(t 3 + 3t 2 u + 4£m 2 + 2m 3 ) + t 4 + 2t 3 u + 5£ 2 m 2 + 4£m 3 + 3m 4 ] 

+ s 3 m[s 4 + 2s 3 (t + m) + 3 s 2 {t + m ) 2 + 2s (t + m ) 3 + ( t 2 + tu + u 2 ) 2 ], 

L = —2 QH 3 u 2 — 2 Q 2 st 2 u[s(t + m) +t 2 + tu + 2m 2 ] — s 2 t[s 2 (t + m ) 2 

+ 2 s(t 3 + 2 t 2 u + 2tu 2 + m 3 ) + t 4 + 2t 3 u + 3t 2 M 2 + 2tu 3 + 2m 4 ]. (A. 6 ) 

e + g -> cc[ 3 Sf } ] + 5-: 

F=yF(e + 9 ^cc[ 3 S< 1) ]+ 9 ), 

T=T(e + g—>cc + 9 ) , 

L = L [e + g ^ cc + g) . (A.7) 

e + g -*• cc [ 3 Pj 8) ] + 5 -: 

967 r 2 e 2 aa 2 

M|s 4 t(s + t) 3 (s + m) 4 (£ + m ) 3 ’ 

T = —8Q 8 t 2 [2s 5 t — 2s 4 m 2 + s 3 t 2 {2t — 3m) — s 2 tu(t 2 + 3£m — 3m 2 ) — stu(2t 3 + £ 2 m — tu 2 
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— m 3 ) + t 2 u 3 {2t + u)] + 2Q 6 t[As 7 (t + u) + 4s 6 (5t 2 + 2 u 2 ) + 4s 5 (5t 3 + 6t 2 u — 2£m 2 
+ 4m 3 ) + 2s 4 (9t 4 - 9 t 3 u + 4t 2 u 2 - 6 tu 3 + 4m 4 ) + s 3 (12t 5 - 2 t 4 u - 29tV + 33tV 

— 12hu 4 + 4m 5 ) — s 2 t(2t 5 + 20t 4 M + 25 t 3 u 2 — t 2 u 3 — 16£m 4 + 4m 5 ) — st 2 M(12t 4 + 26£ 3 m 
+ 2t 2 u 2 + tu 3 + m 4 ) — t 3 u 2 (2t 3 — 6t 2 u — lbtu 2 — 7m 3 )] — 2Q 4 [2s 8 (3t 2 + tu — 2m 2 ) 

+ s 7 (21t 3 + 3t 2 M + 10tM 2 - 8 m 3 ) + s 6 (28t 4 + 15 £ 3 m - 15tV + 22tu 3 - 12m 4 ) + s 5 (27t 5 
+ llt 4 M + 8t\ 2 — 20t 2 M 3 + 26 £m 4 — 8m 5 ) + s 4 (16t 6 + 17 L 5 m — 20t 4 M 2 + 2 t 3 u 3 — 33t 2 u 4 
+ 18tM 5 — 4 m 6 ) + s 3 t(2t 6 + 8t 5 u — 3t 4 u 2 — 19t 3 M 3 — 19£ 2 m 4 — 17£m 5 + 6m 6 ) — s 2 t 2 {2t & 

+ 14t 5 M + 30t 4 M 2 + 33 t 3 M 3 + 46t 2 M 4 + 21 £m 5 + 4m 6 ) — 2st 3 u(2t > + 10£ 4 m + 16t 3 M 2 
+ 20t 2 M 3 + 19 £m 4 + 7m 5 ) — 2 t 6 u 2 (t + m) 2 ] + Q 2 s[s s (7t 2 — 5 tu — 12m 2 ) + s 7 (25t 3 + 3t 2 u 

- 18 tu 2 - 36m 3 ) + s 6 (37t 4 + 25 t 3 u - 12iV - 207u 3 - 60 m 4 ) + s 5 (39t 5 + 39t 4 M + 16tV 

- 18t 2 M 3 - 6tM 4 - 60 m 5 ) + s 4 (29t 6 + 83 t 5 u + 72t 4 u 2 + 88t 3 u 3 + 40t 2 M 4 + 22tu 5 - 36m 6 ) 

+ s 3 (9t 7 + 75t 6 M + 148L 5 m 2 + 176t 4 M 3 + 178£ 3 m 4 + 102tV + 22 tu e - 12m 7 ) + s 2 tu(22t 6 
+ 107t 5 M + 199t 4 M 2 + 2m 3 M 3 + I77t 2 u 4 + 73tM 5 + 9 m 6 ) + st 2 M 2 (17t 5 + 69t 4 M + 107£ 3 m 2 
+ 105t 2 M 3 + 71 £m 4 + 21m 5 ) + 4t 5 M 3 (t + m) 2 ] + s 2 (s + u)\Js 7 u(t + m) + s 6 u(25t 2 + 38 tu 
+ 21m 2 ) + s b (2t 4 + 47t 3 M + 88 t 2 u 2 + 78tu 3 + 35m 4 ) + s 4 (4t 5 + 63t 4 M + 132tV 

+ 156t 2 M 3 + 98tM 4 + 35 m 5 ) + s 3 (2t 6 + 47 £ 5 m + 136t 4 M 2 + 190t 3 M 3 + 156tV + 78 tu 5 
+ 21m 6 ) + s 2 M(13t 6 + 70t 5 M + 136t 4 M 2 + 132t 3 M 3 + 88tV + 38tu 5 + 7m 6 ) + s£m 2 (13£ 5 
+ 47t 4 M + 63 t 3 u 2 + 47t 2 M 3 + 25 £m 4 + 7m 5 ) + 2 t 4 u 3 (t + m) 2 ], 

L = —8 Q 6 t 2 [s 5 t — s 4 u 2 — 2 s 3 t 3 + s 2 t 2 u[t + 3m) + stu{2t 3 + t 2 u — tu 2 — m 3 ) — t 2 u 3 (2t + m)] 

+ 2Q 4 t[2s 7 (t + m) + 4s 6 (2t 2 + m 2 ) — s 5 (t 3 — t 2 u + 2 tu 2 — 8m 3 ) — s 4 (21t 4 + 13 £ 3 m 
+ 18t 2 M 2 — 2 tu 3 — 4m 4 ) — 2s 3 (6f 5 + 10t 4 M + t 3 u 2 + 4£ 2 m 3 — 6tu 4 — m 5 ) + 2 s 2 t(t 5 + 10t 4 M 
+ llt 3 M 2 + 9t 2 u 3 + 6tu 4 + 5m 5 ) + st 2 u(12t 4 + 26 t 3 u + 2t 2 u 2 + tu 3 + m 4 ) + t 3 u 2 (2t 3 

— 6t 2 u — 15 tu 2 — 7m 3 )] — 4 Q 2 [s 8 (t 2 — m 2 ) + s 7 (3t 3 + tu 2 — 2m 3 ) — s e (2t 4 + 2 t 3 u + 2 t 2 u 2 

- tu 3 + 3m 4 ) - s 5 (12t 5 + 16t 4 M + 16t 3 M 2 - t 2 u 3 + tu 4 + 2m 5 ) - s 4 (10t 6 + 28t 5 M + 28t 4 u 2 
+ 13t 3 M 3 - 7t 2 u 4 + 3tM 5 + m 6 ) - s 3 t(t 6 + 9 t 5 u + 18t 4 u 2 + t 3 u 3 - 15t 2 u 4 - lOtu 5 + 2m 6 ) 

+ s 2 t 2 (t 6 + 7 t 5 u + 13t 4 M 2 + 18t 3 M 3 + 35 t 2 M 4 + 24tM 5 + 6 m 6 ) + st 3 u(2t 5 + 10t 4 M + 16t 3 M 2 
+ 20t 2 M 3 + 19tM 4 + 7 m 5 ) + t 6 u 2 (t + m) 2 ] — 2s(s + M)[2s 7 M(t + m) + 2s 6 u(3t 2 + 3 tu 

+ 2m 2 ) + s 5 (5£ 4 + 15t 3 M + 18t 2 M 2 + 14tM 3 + 6 m 4 ) + s 4 (15t 5 + 38t 4 M + 53 t 3 u 2 + 40t 2 M 3 
+ 22tu 4 + 4m 5 ) + s 3 (15t 6 + 52t 5 M + 88t 4 u 2 + 81t 3 M 3 + 47t 2 M 4 + 19tM 5 + 2 m 6 ) + s 2 t(5t 6 
+ 32 t 5 u + 78 t 4 u 2 + 90t 3 u 3 + 68t 2 u 4 + 34tM 5 + 9 m 6 ) + st 2 u(7t 5 + 31t 4 M + A7t 3 u 2 + 39t 2 u 3 
+ 23tM 4 + 7 m 5 ) + 2t 5 u 2 (t + m) 2 ]. (A.8) 
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Figure 1: Representative Feynman diagrams for the partonic subprocesses e + a —> e + 
cc[n ] + a, where a = g,q,q and n = 3 S[ l \ 3 S[ 8 \ 3 -Pj 8) . There are six diagrams of the 

type shown in part (a), two ones of the type shown in part (b), and two ones of the type 
shown in part (c). There are two more diagrams that are obtained from the diagrams of 
the type shown in part (b) by replacing the external gluon lines with quark ones. The CS 
process only proceeds through the diagrams shown in part (a). 
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Figure 2: Schematic representation of e + p 
momentum assignments. 
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Figure 3: The Q 2 distribution of J/if; inclusive production in ep DIS measured with the 


HI detector [12[ is compared with the LO NRQCD (upper histogram) and CSM (lower 


histogram) predictions. 
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Figure 4: The p\ distribution oi J /ip inclusive production in ep DIS measured with the 
HI detector [12[ is compared with the LO NRQCD (upper histogram) and CSM (lower 


histogram) predictions. 
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Figure 5: The z distribution of J/-0 inclusive production in ep DIS measured with the 
HI detector [12[ is compared with the LO NRQCD (upper histogram) and CSM (lower 


histogram) predictions. 
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Figure 6: The y* distribution of J/ip inclusive production in ep DIS measured with the 


HI detector [12[ is compared with the LO NRQCD (upper histogram) and CSM (lower 


histogram) predictions. 
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Figure 7: The W distribution of J/rj) inclusive production in ep DIS measured with the 
HI detector [12[ is compared with the LO NRQCD (upper histogram) and CSM (lower 
histogram) predictions. 
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